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Abstract 

In this paper, we present a three-point without memory iterative method based on Kung and 
Traub’s method for solving non-linear equations in one variable. The proposed method has eighth- 
order convergence and costs only four function evaluations each iteration which supports the Kung- 
Traub conjecture on the optimal order of convergence. Consequently, this method possesses very high 
computational efficiency. We present the construction, the convergence analysis, and the numerical 
implementation of the method. Furthermore, comparisons with some other existing optimal eighth- 
order methods concerning accuracy and basins of attraction for several test problems will be given. 
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1 Introduction 

Solving nonlinear equations is a basic and extremely valuable tool in all fields in science and engineer¬ 
ing. One can distinguish between two general approaches for solving nonlinear equations numerically, 
namely, one-step and multi-step methods. Multi-step methods overcome some computational issues 
encountered with one-step iterative methods. Typically they allow us to achieve a greater accuracy 
with the same number of function evaluations. Important aspects related to these methods are order of 
convergence and optimality. Therefore, it is favorable to attain with fixed number of function evalua¬ 
tions each iteration step a convergence order which is as high as possible. A central role in this context 
plays the unproved conjecture by Kung and Traub It states that an optimal multi-step method 
without memory which uses + 1 evaluations could achieve a convergence order of 2 k . Considering this 
conjecture, many optimal two-step and three-step methods have been presented. 

In the recent years, a large number of multi-step methods for finding simple roots x* of a nonlinear 
equation f(x) = 0 with a scalar function f : D cK aR which is defined on an open interval D (or 
/ : D C C C defined on a region D in the complex plane C) have been developed and analyzed 
for improving the convergence order of classical methods. The Newton-Raphson iteration x n +i := 
x n — P r °bably the most widely used algorithm for finding roots. It is of second order and 
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requires two evaluations for each iteration step, one evaluation of / and one of /'. The Newton- 
Raphson iteration is an example of a one-point iteration, i.e., in each iteration step the evaluations are 
taken at a single point. The basic optimality theorem for one-point iterations (see Traub [231 §5.4] or 
an improved proof in M) shows that an analytic one-point iteration based on k evaluations is of order 
at most k. Thus, the Newton-Raphson iteration is an optimal one-point method with k — 2. 

Some well known two-point methods without memory are described e.g. in Jarratt [12], King m, 
and Ostrowski m- Using inverse interpolation, Kung and Traub m constructed two general optimal 
classes without memory. Since then, there have been many attempts to construct optimal multi-point 
methods, utilizing e.g. weight functions, see in particular Here, we 

will construct a class of eighth-order methods free from second order derivatives with efficiency index 
v^8 ~ 1.68179; recall that the efficiency index of an iterative method of order p requiring k function 
evaluations per iteration step is defined by E(k,p) = ^/p, see [18]. 

The paper is organized as follows: Section [2] is devoted to introduce the ideas for the construction of 
the new optimal class of eighth-order methods based on Kung and Traub’s method by using a Newton- 
step and suitable weight functions. In Section [3] we give the details of the new methods and investigate 
the convergence order; this allows to present a class of optimal three-point methods by using suitable 
weight functions. Particularizing the weight functions we construct a three-parametric family of eighth- 
order optimal iterative root-finding methods. By assigning particular values to these parameters we 
propose two examples for this kind of methods. Numerical performance and comparisons with other 
methods are illustrated in Section 01 In Section [5] we approximate and visualize the basins of attraction 
of the proposed method and compare them with several existing methods, both graphically and by 
mean of some numerical measures. Finally, a conclusion is provided in Sectional 


2 Description of the method 


In this section we construct a new optimal three-point class of iterative methods for solving nonlinear 
equations based on Kung and Traub’s method HU. 

The Kung and Traub’s method is given by 

_ f(Xn) 

Vn •— x n \ •> 

/ (*n) 

< ( 2 - 1 ) 

_ f{Xn)f{yn) fiXn) 
n+1 ' VU (f(Xn)-f(yn)Y'nx n y 

In ([21]) and all forthcoming methods, the iteration rule is used for n — 0,1,... and xo denotes an 
initial approximation of the simple root x*. The convergence order of (12. ip is four with three function 
evaluations in each iteration step. Hence, this method is optimal. We intend to increase the order of 
convergence and extend (12.11) by mean of an additional Newton step 


f(Xn)f(y n ) _ f{Xn) 

(f(x n )-f(y n )) 2 f'(x n y 


x n+ 1 • Zn 


f(*n) 

f'(Zn)' 


( 2 . 2 ) 


Method (|2.2|) uses five function evaluations with convergence order eight. Consequently, this method 
is not optimal. In order to decrease the number of function evaluations, we are going to approximate 
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f'(z n ) by an expression based on f(x n ), f{y n ), f(z n ), and f(x n ), namely 

f(x n ) 


f'(z n ) 


J ( t n , u n )G(s n ) 


with t n := u n := jyj—j, s n jpj- and suitable functions J and G. 

Therefore, we have 

f(x n ) 

y n •— x Ti 


z n : — Vn 


f'{? n y 

f(x n )f(y n ) f(x n ) 


(/(*n) - f(yn)) 2 f'(x n ) 
Xn+1 ' = Zn "777 n n )G(s n ), 

< J \ x n) 


as iteration rule. 


(2.3) 


3 Convergence analysis 

In the following theorem, we analyze the convergence order of method (12.31) . In particular, we find the 
requisites to the weight functions J and G in (|2.3j) guarantee the requested order eight. Although we 
enunciate the method for real functions and a real root, the same can be written (with an identical 
proof) if we have a complex function f : D cC-^C with a complex root x* £ D. 

Theorem 1. Let f : D C M M be an eight times continuously differentiable function with a simple 
zero x* G D, and let J : M 2 M and G:RgK sufficiently differentiable functions in a neighborhood of 
the origin. If the initial point xo is sufficiently close to x*. Then, the method defined by (12.31) converges 
to x* with order eight if the conditions 

Jo,o = Ji,o 2 , t/ 2,0 = 8 , Jo,i = 2 , Js : o — 36, 

and 

Go = 1, G\ = 1, 

with Jij = ^^|( t , u)= (o,o) and G t = ^H| s=0 are fulfilled. 

f(n )( x *\ 

Proof. Let e n := x n — x*, e n ^ y := y n — x*, e UiZ := z n — x* and c n - r J for n G N. Using the fact 
that f(x*) = 0, the Taylor expansion of / at x* yields 

f(x n ) = fix*) (e n + c 2 el + c 3 el H-f c 8 e®) + 0(e®) (3.1) 

and 

f'(x n ) = f'(x*) (l + 2 c 2 e n + 3c 3 e^ + 4c 4 e^ 4 -b 9c 9 e^) + 0(e g n ). ( 3 . 2 ) 

Therefore, we have 

jfj—\ = e n- c 2 e 2 n + (2c| - 2c 3 ) el + (-4c| + 7c 2 c 3 - 3c 4 )e* 

J \ x n) 

+ ( 8 c\ - 20 c| c 3 + 6 c 2 + 10 c 2 c 4 - 4cs)e® 

+ (—16c| + 52c2C 3 — 28 C 2 C 4 + 17c 3 c 4 - c 2 (33c| - 13c 5 ))e® + 0(e 7 n ), 

and 

e n ,y = y n - X* = c 2 e 2 n + (—2c| + 2c 3 )e^ + (4$ - 7c 2 c 3 + 3c 4 )4 
+ (— 8cj + 20c|c 3 — 6 c 2 — 10c 2 c 4 + 4c^)el 

+ (16c| — 52c2C 3 + 28 C 2 C 4 — 17c 3 c 4 + c 2 (33c 3 — 13c5))e^ + O(e^). 
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(3.3) 


We have for f(y n ) also 

f(y n ) = f'(x*) (e n ,y + c 2 e 2 n ^y + c 3 e\ y + • • • + c 8 e^ y ) + 0(e® )W ). 
Therefore, by substituting (13.11) . (13.2j) . and (13.31) into (I2.2jl . we get 

= Z n -x* = (2c\ - c 2 c 3 )e^ - 2(5c| - 7c^c 3 + c\ + c 2 c 4 )e^ 

+ (31c2 — 726 c2C 3 + 2 IC 2 C 4 — 7c 3 c 4 + c 2 (30c 3 — 3c 3 )) + O(e^). 


We get for f(z n ) also 

f( z n) = f\ x *) ( e n,z + c 2 e^ + C 3 e^ H-h C8 e n, z ) + 0(e^ z ). 

From (13.11) and (13.31) . we have 

t n = 77^4 = c 2 e n + (—3c! + 2c 3) e l + ( 8c 2 - 10c 2 c 3 + 3c 4 )e^ 

/ \ x n) 

+ (-20c! + 37 c!c 3 - 8c 3 - 14 c 2C 4 + 4cs)e^ 

+ (48c! - H8c!c 3 + 51c!c 4 - 22c 3 c 4 + c 2 (55c| - 18c 5 )) e b n + 0(e®), 
and from (13.1 1) and (13.41) . we obtain 

u n = 77^4 = (2c| - c 2 c 3 )el + (-12 4 + 15c!c 3 - 2c| - 2c 2 c 4 )e„ 

J \ x n) 

+ (43c! - 89c!c 3 + 23c!c 4 - 7c 3 c 4 + c 2 (33c 3 - 3c 5 ))e® + 0(e®). 


We get from (13.31) and (13.41) 
f(z 1 

Sn = -Tr~\ = ( 2C 2 “ C 3 ) e n “ 2 ( 3c 2 “ 4 C 2 C 3 + C 4 )4 
/ yVn) 

+ (9c^ — 25c! c 3 + 7 c| + llc 2 C 4 — 3c5)e^ 

+ 2 (cl — 18c!c 3 + 15 c!c 4 — 9c 3 c 4 + c 2 (16c| — 7cs)) + 0(e®). 


Expanding J at (0,0) and G at 0 yields 


J{tni u n ) — J 0 ,0 + u nJo,l + + —^^3,0 + 0 (tU ^), 

G(s n ) = Gq + S n G\ + O(s^). 


Substituting (I3.1D - (I3.9D into (12.31) . we obtain 

e n+l = x n+ 1 — x* = i? 4 e n + + Rft&n + ^?7 e n + ^8 e n + O(e^), 


(3-4) 


(3.5) 


(3.6) 


(3.7) 


(3.8) 

(3.9) 


where 

Ra = —(2c! — c 2 c 3 )(—1 + GoJo t o), 

R 5 = —c!(2c! — c 3 )(—2 + Ji ;0 ), 

R& = — 2 C 2( 2c 2 — c 3) (—2c 3 (—1 + Gi) + c!(—12 + 4Gi + J 2j o)) > 

R 7 = — —c!(2c! — c 3 ) (—6c 3 (—2 + J 0 ,i) + c|(—60 + 12Jo,i + ^ 3 , 0 )) • 

By setting i? 4 = • • • = R 7 = 0 and R 8 / 0, the convergence order becomes eight. Obviously, we have 

Jo,o = 1, Go = 1 => R 4 = 0, 

Ji,0 = 2, =>■ = 0, 

^ 2,0 = 8, Gi = 1 =>■ Rq = 0 , 

Jo,i = 2, J 3; o = 36 =>■ R 7 = 0. 
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Consequently, the error equation becomes in this case 

e n+ 1 — ( c 2 ( 2 c^ — C 3 ) (23c| — I 2 C 2 C 3 + C 3 + C 2 C 4 )) + O(e^) 

which finishes the proof of the theorem. □ 


In what follows, we give some concrete explicit representations of (|2.3I) by choosing different weight 
functions satisfying the required conditions for the weight functions J(t n ,u n ) and G(s n ) of Theorem [lj 
We can choose the weight functions J(t n ,u n ) and G(s n ) as 


J(tm Ujj 


1 + at n + (2 + b)u n + (2a + 1 4 - 4a£^ 

1 + (a - 2 )t n + bu n + *2 


(3.10) 


and 


G(s n ) 


1 + cs n 
1 + (c — l)s n 


(3.11) 


with arbitrary a, 6, c G C. It is a simple task to check that the functions J(£ n , u n ) and G(s n ) in (13.lOj) 
and (13. 11 j) satisfy the assumptions of Theorem[1]for all choices of a, 6, c. Hence, three-parametric family 
of optimal eighth-order iterative root-finding methods is obtained. 

By fixing the particular parameters a, 6, c, we are going to give two examples of this family of 
methods. 


Method 1: Set a — b — c— -. Then, we get 


Un •— %r 


Un 


%n+l •— Zr\ 


fjXn) 

f'(Xn) ’ 

f(x n )f(yn 


f (X 


(/(*n) ~ f{Vn )) 2 f{.X n y 

f( z n) (2 + t n + 5?/, n + 4^ + 4/^ 2 + .sv 


f'(x n ) 

v = ±^11 s = 

f{xn) : n /(in ) 1 n /(jfa) 


2 - 3/ n + U n + 2/2 


with t„ = .. - ^ ^ 


Z “h 1 Z — 1 

Method 2: Set a = —-—, b = 1 + z, and c =-. So, we have 


Un •— 


• Vn 


2 

fjXn ) 

f'(x n ) ’ 

f(x n )f{y r 


f (Xr 


%n -\-1 •— ^77 


(/fan) - /fan)) 2 /'fan)’ 

/fan) ( \ + (^ 2 ^)^n + (i + 3)tt n + (i + 2)/% + 4(^jp)/3 1 + (^far; 

/'fan) l 


1 + (^ 2 “)t n + (z + 1)^71 + 


i + (V) s . 


(3.12) 


(3.13) 


with / = v = lilzil s = /( z ") 

n f(x n )> n f(x n )’ n f(yn) m 

We will apply in the next sections the new methods (13.121) and (I3.13j) to several benchmark examples 
and will compare the new methods with some existing optimal three-point methods of order eight having 
the same optimal computational efficiency index equal to y/8 1.68179, see Hi. 

The existing methods that we are going to use to compare are the following: 
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Method 3: The method by Chun and Lee [5] is given by 


Un •— %n 


z n : — Un 


x n +1 : 


— Zn. 


fjXn) 

f'(Xn ) ’ 

fiVn) _ 

/'On) 

/On) _1_ 

/'On) ' (1 - H(t n ) ~ J(s n ) - P( % )) 2 


(1 _ /(y«) \ 

V 1 f(xn)J 


(3.14) 


with weight functions 


t 2 t 3 s u 

ff(*n) = -/3-7 + *n + f-f, J(*n) = /3+f, %) = 7+y, 

where t n — Sn = j^y, u n = , and / 3,7 E M. Note that the parameters (3 and 7 cancel when 

used in f)3.14p . Hence, their choice has no contribution to the method. 


Method 4: The method by B. Neta on, see also formula (9)], is given by 

_ _ /On) 

Vn '- Xn r(x n y 

< „ _ „ /On) + Af(y n ) f(y n ) A ^ (3.15) 

n - Jn /(x„) + (A-2)/( ita )'/'(x„)’ ’ 

k z n+ i := y n + Sif 2 (x n ) + 5 2 f 3 (x n ), 

where 


^y = f(Vn) -/On), 


Cy ~ ^ 


1 (Vn %r 


K 

62 = - 


y \ ^y 
Cy ~ Cz 
r y r z 


/'On )J ’ 


F z = f(z n ) - /On), 

1 2^77 3^71 

Q = tt - 


F z V F* 


/'(*«)/ ’ 


0 — Cy + S 2 Fy. 


We will use A = 0 in the numerical experiments of this paper. 


Method 5: The Sharma and Sharma method [22j is given by 

/On) 


Vn := X r 


/'On) ’ 


/On 


z :=y _1M_ _ 

f f (x n ) f(x n ) - 2f{y n y 


%n -\-1 •— 


/Qn^nl/Qn) 

f[Xn,Z n }f[y n ,Z n } 


W(t n ), 


with the weight function 


W(t n ) = 1 + 


1 T Qtti 


-, o E M., 


and t n = jTj- We will use a = 1 in the numerical experiments of this paper. 


(3.16) 
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Method 6: 


The method from Babajee, Cordero, Soleymani and Torregrosa [2] is given by 


Un •— x n 


f( x n) (■,,( /W yt 

f'M \ \f'M) J ’ 


< 


Zn •— Vn 


f(Vn) 



f{yn) \ 

f( x n)J 


, ( /On) A 2 - ( f(yn ) \ 4 f(Zn) 

/(^n) _ \f( x n)J \f( x n)J f(y n ) 

L fiVn) f(z n ) \ 2 

V f( x n) f( x n)J 


(3.17) 


4 Numerical examples 

The particular cases (13.121) and (13.131) of the the three-point method (12.3j) are tested on a number 
of nonlinear equations. To obtain a high accuracy and avoid the loss of significant digits, we em¬ 
ployed multi-precision arithmetic with 20 000 significant decimal digits in the programming package 
Mathematica. 

In order to test our proposed methods (13.121) and (13.131) . and also to compare them with the methods 
(I3.14|) . (13.151) . (I3.16|) . and (13.171) . we compute the error, the computational order of convergence (COC) 
by the approximate formula m 


In K^+i -x*)/(x n -x*)| 
In | {x n -x*)/(x n -i -x *)\’ 


and the approximated computational order of convergence (ACOC) by the formula [H] 


(4.1) 


ACOC « 


In \(x n +l /\%n %n— 1 )| 

In \{x n - x n -i)/{x n -i - x n - 2 )[ 


(4.2) 


It is worth noting that COC has been used in the recent years. Nevertheless, ACOC is more 
practical because it does not require to know the root x*. For a comparison among several convergence 
orders, see [H]. Moreover, we should note that the results for these formulas not always coincide with or 
approximate the exact convergence order of the method when they are applied to a particular example. 
The reason is that we have in the error equations of the methods some coefficients that depend on c& 
(see the proof of Theorem [Tj) . Hence, these c/ds may vanish or vary for different kinds of examples. 
But, in general, a “random” example should provide good approximations for the order of convergence 
of the method. 

On the other hand, it is nice to note that, given an iterative method, computing COC or ACOC on 
several examples is a good experiment to check theoretical errors in the deduction of the method and 
to check practical errors in the implementation of the method in a computer. For a general problem it 
will be difficult that COC or ACOC approach the theoretical order of convergence by chance. 

We have used both COC and ACOC for checking the accuracy of the considered methods. Note 
that both COC and ACOC give already for small values of n good experimental approximations to 
convergence order. 

In what follows, we are going to perform this kind of numerical experiments with the four test 
functions fj(x ), j = 1,...,4, that appear in Table HJ In every case, and using the six eighth-order 
iterative methods described in the paper, we are going to reach the root x* starting in the point xq. 

In Table O our new three-point methods (13.121) and (13.131) are tested on the four nonlinear equations 
fj(x ) = 0, j = 1,2, 3,4, and compared them with the methods (13.141) . (13.151) . (13.161) . and (13.171) from 
other authors. We abbreviate (13.12l) - (13.17l) as M1-M6. Notice that, to estimate the COC and the 
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test function fj 

root x* 

initial guess xq 

fi(x) = ln(l + x 2 ) + e x 2 ~3x s j nx 

0 

0.35 

f 2 (x) = 1 + e 2+x ~ x<2 + x 3 — cos(l + x) 

-1 

-0.3 

h(x) = (1 + z 2 )cos " + ln( ’i+5 +2) 

-1 

-1.1 

fa{x) = x 4 + sin y — 5 

V2 

1.5 

Table 1: Test functions /i,..., /p root x*, and initial guess xq. 


ACOC, it has been enough to use n — 3 in (14.11) and (14.21) to get excellent approximations of the order 
of convergence. 



Ml 

M2 

M3 

M4 

M5 

M6 

fi, x 0 = 0.35 

|a?i - x* 

0.140e—3 

0.318e—3 

0.721e—4 

0.893e—4 

0.753e—4 

0.347e—3 

\X2 ~ X*\ 

0.583e—28 

0.562e—25 

0.230e—30 

0.126e—30 

0.619e—31 

0.471e—25 

\x 3 - X*\ 

0.362e—223 

0.531e—199 

0.252e—242 

0.200e—245 

0.128e—247 

0.546e—200 

COC 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

ACOC 

7.9999 

8.0000 

7.9999 

7.9999 

7.9999 

7.9999 

h, xo = -0.3 

|^i — x*\ 

0.526e—4 

0.113e—3 

0.157e—3 

0.763e—4 

0.871e—4 

0.411e—3 

\X2 ~ X*\ 

0.534e—37 

0.263e—33 

0.119e—33 

0.540e—35 

0.134e—34 

0.377e—29 

a; 3 - x* \ 

0.599e—301 

0.226e—270 

0.138e—274 

0.342e—284 

0.438e—281 

0.189e—237 

COC 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

ACOC 

7.9999 

7.9999 

7.9998 

7.9999 

7.9999 

7.9999 

/3, X 0 = -1.1 

\xi - x*\ 

0.235e—7 

0.298e—7 

0.614e—8 

0.388e—8 

0.175e—8 

0.554e—8 

\x 2 - X*\ 

0.393e—60 

0.373e—59 

0.328e—65 

0.254e—67 

0.154e—70 

0.426e-66 

a; 3 - x*\ 

0.239e—482 

0.222e—474 

0.217e—523 

0.877e—541 

0.5821e—567 

0.528e—531 

COC 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

ACOC 

7.9999 

7.9999 

8.0000 

7.9999 

8.0000 

8.0000 

fa, Xq = 1-5 

\xi - x*\ 

0.286e—8 

0.602e—8 

0.433e—8 

0.327e—10 

0.642e—10 

0.281e—8 

\x 2 - x*\ 

0.108e—68 

0.181e—65 

0.134e—66 

0.369e—84 

O.lOle—81 

0.341e—68 

\x 3 - x*\ 

0.460e—552 

0.121e—525 

0.116e—534 

0.967e—676 

0.389e—656 

0.161e—547 

COC 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

8.0000 

ACOC 

8.0000 

8.0000 

7.9999 

7.9999 

7.9999 

8.0000 


Table 2: Errors, COC, and ACOC for the iterative methods (13.12l) - (13.17l) (abbreviated as M1-M6) 
applied to the find the root of test functions fi ,... , given in Table [lj 


5 Dynamic behavior 

We already observed that all methods converge if the initial guess is chosen suitably. We now investigate 
the regions where we must choose the initial point to achieve the root. In other words, we numerically 
approximate the domain of attraction of the zeros as a qualitative measure of how demanding is the 
method on the initial approximation of the root. To answer the important question on the dynamical 
behavior of the algorithms, we investigate the dynamics of the new methods (13.121) and (13.131) and 
compare with common and well-performing methods from the literature, namely (13.141) . (13.151) . (13.161) . 

























and (I3.17p . We recall in the following line some basic concepts such as basin of attraction. For more 
details and many other examples of the study of the dynamic behavior for iterative methods, one can 

consult Pfl El SI El El Um mi 1231 E5J. 

Let Q : C —>> C be a rational map on the complex plane. For z G C, we define its orbit 
as the set orb(z) = {z, Q(z ), Q 2 (z), ... }. A point zq G C is called periodic point with minimal 
period m if Q m (zo ) = z$ where m is the smallest positive integer with this property (and thus 
{^o, Q(^o), • • •, Q m_1 (^o)} is a cycle). A periodic point with minimal period 1 is called fixed point. 
Moreover, a fixed point z$ is called attracting if \Q'(zo)\ < 1, repelling if \Q'(zo)\ > 1, and neutral 
otherwise. The Julia set of a nonlinear map Q(z) : denoted by J(Q), is the closure of the set of its 
repelling periodic points. The complement of J(Q ) is the Fatou set F(Q). 

In our case, the six methods fl3.12l) - fl3.17l) provide iterative rational maps Q(z) when they are 
applied to find the roots of complex polynomials p(z). In particular, we are interesting in the basins of 
attraction of the roots of the polynomials where the basin of attraction of a root z* is the complex set 
{^o G C : orb(zo) —^ z*}. It is well known that the basins of attraction of the different roots lie in the 
Fatou set F(Q). The Julia set J(Q) is, in general, a fractal and, in it, the rational map Q is unstable. 

For the dynamical and graphical point of view, we take a 512 x 512 grid of the square [—3, 3] x 
[—3,3] C C and assign a color to each point zo G D according to the simple root to which the 
corresponding orbit of the iterative method starting from zq converges, and we mark the point as black 
if the orbit does not converge to a root in the sense that after at most 15 iterations it has a distance 
to any of the roots which is larger than 10 -3 . We have used only 15 iterations because we are using 
eighth-order methods so, if the method converges, it is usually very fast. In this way, we distinguish 
the attraction basins by their color. 


Test polynomials 

Pl(z) = z 2 - 1 

p 2 (z) =z 3 -z 

p 3 (z) = z(z 2 + l)(z 2 + 4) 

p 4 (z) = (z 4 - 1 )0 2 + 2 i) 

p 5 (z) = z 7 - 1 

p 6 (z) = (10z 5 - 1 )(z 5 + 10) 


Roots 

1, -1 
0 , 1 , -1 

0, 2 i, —2 i, i, —i 


1, i, —1, —i, —1 + i, 1—i 

e 2kni / 7 , k = 0,..., 6 

(l^) 1/5 e 2W5 , (—lO) 1 /^/ 5 , k = 0,..., 4 


Table 3: Test polynomials pi(z), ... ,ps(z) and their roots. 

Different colors are used for different roots. In the basins of attraction, the number of iterations 
needed to achieve the root is shown by the brightness. Brighter color means less iteration steps. Note 
that black color denotes lack of convergence to any of the roots. This happens, in particular, when 
the method converges to a fixed point that is not a root or if it ends in a periodic cycle or at infinity. 
Actually and although we have not done it in this paper, infinity can be considered an ordinary point 
if we consider the Riemann sphere instead of the complex plane. In this case, we can assign a new 
“ordinary color” for the basin of attraction of infinity. Details for this idea can be found in mi- 

Basins of attraction for the six methods (I3.12p - fl3.17p for the six test problems Pi(z) = 0, i = 1,..., 6, 
are illustrated in Figures [l]-[6] from left to right and from top to bottom. 

From the pictures, we can easily judge the behavior and suitability of any method depending on 
the circumstances. If we choose an initial point zo in a zone where different basins of attraction touch 
each other, it is impossible to predict which root is going to be reached by the iterative method that 
starts in z$. Hence, it is not a good choice. Both the black zones and the zones with a lot of colors are 
not suitable to take the initial guess zq when we want to achieve a precise root. The most attractive 
pictures appear when we have very intricate frontiers between basins of attraction and they correspond 
to the cases where the method is more demanding with respect to the initial point and its dynamic 
behavior is more unpredictable. 
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Figure 3: Comparison of basins of attraction of methods (I3.12l) - (j3.17l) for the test problem 7 ) 3 ( 2 ) = 
z(z 2 + 1 )( 2 2 + 4) = 0. 



Figure 4: Comparison of basins of attraction of methods (I3.12l) - (j3.17l) for the test problem 7 ) 4 ( 2 ) = 
(z 4 - l)(z 2 + 2i) = 0 . 
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Figure 5: Comparison of basins of attraction of methods (j3.12|) - (|3.17|) for the test problem p$(z) = 
z 7 - 1 = 0. 



Figure 6: Comparison of basins of attraction of methods (I3.12IMI3.17D for the test problem po(z) = 
(10z 5 - 1 )(z 5 + 10) = 0. 


12 













































































































Figure 7: Basins of attraction of the method of Theorem[l]to solve po(z) = (10z 5 — 1 )(z 5 +10) = 0 with 
different choices of the parameters a, b and c. Left: a — b — c— —1. Center: a — b— —1, c = —1 + i. 
Right: a —1/2, b = —1, c = — 2 + i. 

The first and the second graphics in the Figures [THB] correspond to the same general method (|2.3j) 
with different choices of the parameters a, 6 , c in (13.101) and (13. lip . We can see from the graphics in 
Figure Hand the first two pictures of Figure [ 6 ] that even small changes in the parameters may lead to 
completely different behaviors. 

Finally, we have included in Table 0] the results of some numerical experiments to measure the 
behavior of the six iterative methods (13.12l) - (13.17l) in finding the roots of the test polynomials Pj(z), 
j — 1,..., 6 . To compute the data of this table, we have applied the six methods to the six polynomials, 
starting at an initial points zo on a 512 x 512 grid in the rectangle [—3, 3] x [—3, 3] of the complex plane. 
The same way was used in Figures [l]-[ 6 ] to show the basins of attraction of the roots. In particular, 
we decide again that an initial point z$ has reached a root z* when its distance to z* is less than 
10 -3 (in this case zo is in the basin of attraction of z*) and we decide that the method starting in zo 
diverges when no root is found in a maximum of 15 iterations of the method (in this case, we say that 
zo is a “nonconvergent point”). In Tabled! we have abbreviated the methods (I3.12I) - (I3.17I) as M1-M6, 
respectively. The column I/P shows the mean of iterations per point until the algorithm decides that a 
root has been reached or the point is declared nonconvergent. The column NC shows the percentage of 
nonconvergent points, indicated as black zones in the figures. It is clear that the nonconvergent points 
have a great influence on the values of I/P since these points contribute always with the maximum 
number of 15 allowed iterations, In contrast, “convergent points” are reached usually very fast due to 
the fact that we are dealing with eighth-order methods. To reduce the effect of nonconvergent points, 
we have included the column Iq/C which shows the mean number of iterations per convergent point. If 
we use either the columns I/P or the column Iq/C to compare the performance of the iterative methods, 
we clearly obtain different conclusions. 

6 Conclusion 

We have introduced a new optimal class of three-point methods without memory for approximating a 
simple root of a given nonlinear equation which use only four function evaluations each iteration and 
result in a method of convergence order eight. Therefore, the Kung and Traub’s conjecture is supported. 
Numerical examples and comparisons with some existing eighth-order methods are included and confirm 
the theoretical results. The numerical experience suggests that the new class is a valuable alternative 
for solving these problems and finding simple roots. We used the basins of attraction for comparing 
the iterative algorithms and we have included some tables with comparative results. 


13 



























Polynomial 

Method 

I/P 

NC (%) 

Ic/C 

Polynomial 

Method 

I/P 

NC (%) 

Ic/C 

M z ) 

Ml 

2.53 

0.244 

2.50 

Pi(z) 

Ml 

6.85 

24.7 

4.17 


M2 

2.29 

0.00798 

2.28 


M2 

6.48 

22.0 

4.07 


M3 

2.20 

0.195 

2.18 


M3 

4.07 

0.888 

3.97 


M4 

2.17 

0.195 

2.15 


M4 

4.21 

1.84 

4.01 


M5 

2.13 

0.195 

2.10 


M5 

3.95 

4.40 

3.44 


M 6 

6.01 

70.9 

2.09 


M 6 

4.45 

20.1 

3.56 

P2(z) 

Ml 

3.54 

0.798 

3.45 

Ph(z) 

Ml 

7.27 

27.0 

4.42 


M2 

3.10 

0.340 

3.06 


M2 

7.00 

25.2 

4.30 


M3 

2.88 

0 . 

2.88 


M3 

4.81 

3.36 

4.45 


M4 

2.82 

0.00456 

2.82 


M4 

5.07 

5.71 

4.47 


M5 

2.73 

0 . 

2.73 


M5 

4.59 

7.04 

3.80 


M 6 

4.32 

27.6 

2.81 


M 6 

5.03 

21.4 

4.02 

Pz{z) 

Ml 

3.88 

3.57 

3.47 

p&{z) 

Ml 

7.36 

24.4 

4.90 


M2 

3.57 

2.19 

3.31 


M2 

6.96 

21.7 

4.73 


M3 

2.99 

0.0122 

2.99 


M3 

4.69 

2.33 

4.44 


M4 

2.94 

0.0334 

2.94 


M4 

4.89 

4.03 

4.46 


M5 

2.82 

0 . 

2.82 


M5 

4.44 

3.98 

4.01 


M 6 

3.28 

5.46 

2.99 


M 6 

5.26 

11.9 

4.70 


Table 4: Measures of convergence of the iterative methods (I3.12p ~ fl3.17l) (abbreviated as M1-M6) 
applied to find the roots of the polynomials Pj(z), j = 1,..., 6 . 
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